Surface Enhanced Raman Spectroscopy of Graphene 
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Surface enhanced Raman scattering (SERS) exploits surface plasmons induced by the incident 
field in metallic nanostructures to significantly increase the Raman intensity. Graphene provides 
the ideal prototype two dimensional (2d) test material to investigate SERS. Its Raman spectrum 
is well known, graphene samples are entirely reproducible, height controllable down to the atomic 
scale, and can be made virtually defect-free. We report SERS from graphene, by depositing arrays 
of Au particles of well defined dimensions on graphene/SiO2(300nm)/Si. We detect significant en- 
hancements at 633nm. To elucidate the physics of SERS, we develop a quantitative analytical and 
numerical theory. The 2d nature of graphene allows for a closed-form description of the Raman 
enhancement. This scales with the nanoparticle cross section, the fourth power of the Mie enhance- 
ment, and is inversely proportional to the tenth power of the separation between graphene and the 
nanoparticle. One consequence is that metallic nanodisks are an ideal embodiment for SERS in 2d. 
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INTRODUCTION 

Graphene is at the center of a significant research 
effort [1-6]. Near-ballistic transport at room temper- 
ature and high mobility [5-11] make it interesting for 
nanoelectronics [12-15], especially for high frequency 
applications[16]. Furthermore, its transparency and me- 
chanical properties are ideal for micro and nanomechan- 
ical systems, thin-film transistors, transparent and con- 
ductive electrodes, and photonics [17-24]. Graphene is 
also an ideal test-bed for some long-standing problems, 
such as the Raman spectra of carbon materials. Here we 
show that this conceptually simple material (due to its 
low-dimensionality) can be helpful in understanding the 
basics of Surface Enhanced Raman Scattering (SERS). 

Graphene layers can be identified by inelastic [25] and 
elastic light scattering [26, 27]. Raman spectroscopy 
allows monitoring of doping, defects, strain, disorder, 
chemical modifications and edges[25, 28-41]. The Ra- 
man signal can be enhanced for flakes deposited on cer- 
tain substrates, such as the common Si-|-Si02, due to 
interference in the Si02 layer, resulting into enhanced 
field amplitudes within graphene [42-45]. 

Another way to increase the Raman signal is to 
perform SERS experiments[46, 47]. SERS is widely 
used [48-50], and enhancements as large as 14 orders of 
magnitude can be achieved (enough for single-molecule 
detection[51]). However, even though the technique is 
more than 30 years old [50], the exact nature of SERS 
is still debated[48]. Furthermore, the particular mecha- 
nism might be different, depending on whether the Ra- 
man processes involved are resonant or not. In princi- 
ple, even a single metallic nanostructure, e.g., a metallic 
nanotip, can induce SERS at its apex, giving rise to the 
so-called tip-enhanced Raman scattering (TERS) [52-54]. 
The most important feature that makes TERS so attrac- 
tive is its capability of optical sensing with high spatial 



resolution beyond the light diffraction limits[55, 56]. 

Most SERS-active systems studied to-date are based 
on random nanostructures, whose properties vary from 
experiment to experiment making quantitative compari- 
son between theory and experiment difficult. Graphene 
offers a unique model system where SERS effects could be 
studied in detail. Its Raman spectrum is well known, be- 
ing investigated in several hundreds papers in the past 4 
years. Graphene samples are very reproducible and offer 
an atomic-precision control on the number of layers, thus 
allowing a smooth transition from a purely 2d case to a 
3d one. Furthermore, as both resonant and non-resonant 
Raman scattering can be in principle possible, such as in 
chemically modified graphene[33], a distinction between 
different enhancement mechanisms could be made. Here 
we focus on the resonant case, where we believe the en- 
hancement is mostly due to near-field plasmonic effects 
in the vicinity of metal particles [46, 52]. 



EXPERIMENTAL 

Graphene flakes are prepared on Si+300nni Si02 
by micromechanical cleavage [1]. Single layer graphene 
(SLG) is identified by a combination of optical 
contrast [26, 27] and Raman spectroscopy [25]. Electron 
beam lithography in combination with thin metallic film 
deposition (5nm Cr-|- 80nm Au) and lift-off are utilized 
to prepare three sets of metallic dots, as well as a set of 
contacts for transport measurements. Fig 1. One set is 
placed directly on top of graphene, one partially covers 
graphene and partially rests on Si02, and the last com- 
pletely on Si02. The dots sizes and the configurations of 
the arrays can be seen on Fig.l. During lift-off, metallic 
dots are slightly shifted from their lithographically de- 
fined positions, probably by capillary forces. This indi- 
cates poor adhesion of Cr/Au dots onto graphene. Note 
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FIG. 1: Scanning Electron Microscopy images (in false 
colours) of our SERS sample. False colors denote: purple 
- Si02; bluish - graphene; yellow - Au electrodes and dots, a) 
An overall image of the sample. b,c) Golden dots on Si02 



that the dots on Si02 still occupy the positions defined 
by the lithography procedure. 

Raman spectra are recorded with a Renishaw RMIOOO 
spectrometer, equipped with a piezoelectric stage (PI) 
able to shift the sample at nanometer steps. Line scans 
are recorded across the patterned arrays, as shown in 
Figs 2,3, for 488, 514 and 633nm excitation. 



RESULTS AND DISCUSSION 

The Raman enhancement is defined as the ratio of the 
Raman intensity measured on the graphene covered by 
dots, compared to that measured outside the dots, but 
still on the graphene layer. Fig. 2. 

Fig. 3 shows representative Raman spectra measured 
at 633nm. A clear enhancement is seen when comparing 
the patterned and the unpatterned graphene. 

The Raman spectrum of graphene consists of a set of 
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FIG. 2: Raman linescans. Dotted red lines indicate the 
graphene edge. The intensity ratios are a) large dots: Point 
1/Point 3; Point 2/Point 3; b) small dots: Point 1/Point 3; 
Point 2/Point 3; c) large dots: Point 4/ Point 6; Point 5/Point 
6; small dots: Point 3/Point 1; Point 2/Point 1 



distinct peaks. The G and D appear around 1580 and 
1350 cm^^, respectively. The G peak corresponds to the 
E2g phonon at the Brillouin zone center (F point). The D 
peak is due to the breathing modes of six-atom rings and 
requires a defect for its activation [40, 57, 58]. It comes 
from TO phonons around the K point [40, 58], is active 
by double resonance (DR)[57], and is strongly dispersive 
with excitation energy due to a Kohn Anomaly at K[31]. 
DR can also happen as intra-valley process, i. e. con- 
necting two points belonging to the same cone around 
K (or K'). This gives the so-called D'peak, which is 
at~ 1620 cm^^ in defected graphite measured at 514nm. 
The 2D peak is the second order of the D peak. This is 
a single peak in single layer graphene (SLG), whereas it 
splits in four in bilayer graphene (BLG), reflecting the 
evolution of the band structure[25]. The 2D' peak is the 
second order of D'. Since 2D and 2D' originate from a 
process where momentum conservation is satisfied by two 
phonons with opposite wavevectors, no defects are re- 
quired for their activation, and are thus always present. 
Each Raman peak is characterized by its position, width. 
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FIG. 3: Representative Raman spectra measured across a line 
scan moving from outside to inside the Au dot patterned area 
for 633 excitation. A clear enhancement of all peaks is seen 
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FIG. 4: Ratio of height and areas of the G and 2D peaks mea- 
sured on the patterned regions compared to those measured 
outside, as a function of the excitation energy 



height, and area. The frequency-integrated area under 
each peak represents the probabihty of the whole process. 
It is more robust with respect to various perturbations 
of the phonon states than width and height [59]. Indeed, 
for an ideal case of dispersionless undamped phonons 
with frequency Wph, the shape of the n-phonon peak is a 
Dirac S distribution oc S{ijj — nojph), with zero width, in- 



finite height, but well-defined area. If the phonons decay 
(e. g, into other phonons, due to anharmonicity, or into 
electron- hole pairs, due to electron-phonon coupling), the 
6 lineshape broadens into a Lorentzian, but the area is 
preserved, as the total number of phonon states cannot be 
changed by such perturbations. If phonons have a weak 
dispersion, states with different momenta contribute at 



slightly different frequencies. This may result in an over- 
all shift and a non-trivial peak shape, but frequency inte- 
gration across the peak means counting all phonon states, 
as in the dispersionless case. Thus, the peak area is pre- 
served, as long as the Raman matrix element itself is 
not changed significantly by the perturbation. The lat- 
ter holds when the perturbation (phonon broadening or 
dispersion) is smaller than the typical energy scale deter- 
mining the matrix element. Converting this into a time 
scale using the uncertainty principle we have that, if the 
Raman process is faster than the phonon decay, the total 
number of photons emitted within a given peak (i. e., in- 
tegrated over frequency across the peak), is not affected 
by phonon decay, although their spectral distribution can 
be. Even if the graphene phonons giving rise to the D 
and D' peaks are dispersive [31], their relative change with 
respect to the average phonon energy is at most a few %, 
thus we are in the weakly dispersive case. The phonon 
decay in graphene is in the ps timescale, while the Ra- 
man process is faster, in the fs timescale[30, 60, 61]. We 
thus consider both the area, A(2D)/A(G), and height, 
I(2D)/I(G), ratios [59]. Fig 4 plots them for the two dot 
sizes and as a function of excitation energy. 

We model our experiment with the calculation box 
shown in Fig. 5. Starting from the bottom, this con- 
sists of a semi-infinite Si substrate, a SOOnm Si02 and 
SLG of effective thickness 0.335nni. On SLG, we have a 
Au/Cr disk with thickness 80nm/5nm, with the diameter 
set to either 140 or 210nm, according to the experimental 
Au dot size. We time-integrate Maxwell's equations us- 
ing the finite-difference time-domain method(FDTD)[62] 
as implemented in Refs[63, 64] (see Methods). For the 
absorbing boundary conditions in the vertical direction 
we use the perfectly-matched-layer method (PML)[65], 
while in the lateral directions we use periodic boundary 
conditions simulating an infinite two-dimensional square 
array of Au/Cr nanodisks. We previously investigated 
the plasmonic resonances of similar nanoparticles (pre- 
pared in exactly the same conditions as in the present 
work) [66, 67]. This allowed us to extract the Au and Cr 
optical constants as obtained in our evaporators [68]. We 
also recently measured the optical constants of graphene 
by spectroscopic ellipsometry[69]. The dispersive mate- 
rials Au, Cr and graphene are described here by Drude- 
Lorentz models, each fitted to our experimental data. 
These are shown in Fig. 11 in Methods. Finally, for sim- 
plicity we use n=1.46 for Si02, n=4 for Si[70]. 

We only consider normal incidence and emission, rel- 
evant for our Raman backscattering experiments. The 
incident field is a wide spectrum plane wave (i.e. a nar- 
row Gaussian temporal profile) coming from the top. We 
monitor the electric fields E^ (r, t) and Ey (r, t) at each 
grid point on the graphene plane, and upon Fourier trans- 
form we get the tangential-field amplitude £'|| (r, uj) in fre- 
quency domain. This is enhanced compared to the inci- 
dent field, due to substrate interference, and the Surface 
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FIG. 5: Simulation box: 80nm/5nm Au/Cr nanodisks on 
SLG sitting on a Si02/Si. The lateral periodicity is 320nm in 
both directions. We consider nanodisk diameters of 210 and 
140nm, corresponding to the large and small dots 



Plasmon Resonance (SPR) near field of the nanodisks. 

Assuming the graphene absorption at a particular 
point to be proportional to the incident tangential field 
intensity at that point, the Raman emission from a par- 
ticular point to be proportional to the corresponding 
Stokes-shifted intensity, and the emission from points un- 
derneath the nanodisks to be reabsorbed and lost, we 
may approximate the total Raman signal as [46]: 



IsERS oc 



J \E\\ir,,,)\^\E\\ir,iOs)\'dS' (1) 



where the integration is performed over the area not di- 
rectly underneath the nanodisks, and w^ = w — 6oj is the 
Stokes shifted frequency. Slo = 27rcj^, with i^ is the Ra- 
man shift (in cm~^). The outcome of Eq. 1 is normalized 
by the corresponding calculation for suspended graphene 
(where the integral is over all the area since there are 
no nondisks to cover the emission). We note that, given 
the large absorption of ~2.3% per graphene layer [71], this 
approximation becomes questionable for local absorption 
enhancements greater than'^43. To amend this we cut 
above 43. Emission saturation, on the other hand, cannot 
be reached because the Raman efficiency, even though 
larger than other materials due to the process being al- 
ways resonant, is in absolute terms very small [45], and 
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FIG. 6: The total (patterned) enhancement factors for a) the 
G and b) 2D peaks. The dotted hne is the corresponding 
interference (unpatterned) enhancement factor. 



would require ~ 10^^ enhancement to exceed 100%. 

In order to account for the graphene layer thickness 
and still keep the simulation reasonably fast, we employ 
an anisotropic grid with 0.335nm spacing in the vertical 
and 2nm in the lateral dimension. By calculations on 
smaller cells we verified that such a grid introduces small 
errors, typically less than 5% and never exceeding 10%. 
Considering our approximations up to this point (normal 
incidence and the approximation in Eq. 1), we find this 
to be well within the overall simulation errors. 

The 300nm Si02/Si substrate interferometrically in- 
creases not only the visibility[26, 27] but also the Ra- 
man signal of graphene [42, 44, 72]. Here we expect an 
additional enhancement due to the Au nanodisks SPR 
near field. In order to distinguish between the two ef- 
fects (substrate interference and SPR) , we separately cal- 
culate both cases for unpatterned SLG on Si02/Si and 
SLG patterned with the Au nanodisks, still on Si02/Si. 
We define the interference enhancement factor F as the 
ratio of the Raman signal from unpatterned SLG on 
Si02/Si to that of suspended SLG, F = lunpatt/Isusp, 
and the total enhancement factor F' as the ratio of the 
Raman signal from patterned SLG on Si02/Si to that 
of suspended SLG: F' = Ipatt/Isusp = Isers/Iq- Fig. 6 
plots such factors for the G and 2D peaks. This gives a 
maximum interference-related enhancement F at SSOnm 
of~2.5 for the G peak and ~2 for the 2D. However, this 
is very modest compared to the total enhancement F' 
when the nanodisks are taken into account. The total 
enhancement reaches up to 50, and is maximum at dif- 
ferent wavelengths depending on the nanodisk diameter. 
The shoulder at 550nm is likely related to interference, 
because (i) it is at the same frequency as the interference 
peak and (ii) it is stronger the further the plasmon peak. 

There is a different enhancement for the two disk di- 




FIG. 7: Tangential field intensity distribution at 633 nm for 
the a) 210 nm and b) 140 nm nanodisk diameters. 
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FIG. 8: Normalized enhancement factors F' /F for a) the G- 
band and b) the 2D-band. Squares and crosses are the exper- 
imental Raman intensities and areas. 



ameters, not only in peak value, but also in wavelength. 
This is a size effect: the nanodisks are comparable to the 
incident laser wavelength. Thus, retardation effects and 
higher order multipole terms become important, mod- 
ifying the plasmon response as a function of disk size. 
Fig. 7 plots the distribution of the tangential intensity 
enhancement in the graphene plane at 633nm, which is 
approximately proportional to the absorption enhance- 
ment at that wavelength. Different patterns appear for 
the two sizes, demonstrating the point made above. 

Fig. 8 compares F' jF to the corresponding experimen- 
tal Raman intensity ratios. Considering all approxima- 
tions made, and the fact that some of the measurements 
where on top of distorted parts of the nanodisk arrays, 
the agreement is good. We find low enhancement for 488 
and 514nm, and large for 633nm. Also, we reproduce the 
higher enhancement for the G peak in the small disks 
at 633nm, and for the 2D peak in the large disks. The 
quantitative agreement for the 2D peak is not as good, 
however, this is expected given that its intensity signif- 
icantly depends on electron-electron interactions, which 
could change in presence of gold [59]. 

The agreement between experiment and simulation 



is encouraging. It also allows to get new physical in- 
sights into the SERS process of an extended 2d system, 
like graphene. The basic physics and detailed theory of 
the electromagnetic contribution to SERS on adsorbed 
molecules is well known [46]: both absorption and emis- 
sion are enhanced due to interaction with surface plas- 
mons, with an expected overall dependence on the fourth 
power of the SPR-mediated field enhancement [46, 73]. 
However, a detailed theory for 2d systems is lacking. 
Such formulation is challenging for our experiment: the 
Au particles do not have a regular spherical or ellipsoidal 
shape, they are large and so multipoles higher than dipole 
contribute, there is a thin Cr layer, we are on Si02/Si 
giving additional interference and enhancement effects. 
We thus consider the simplified case of regular-shaped 
small Au nanoparticles inside a uniform medium. This 
will provide all the new physics for SERS in graphene, or 
generally any 2d system. The final connection between 
experiment and theory still relies on simulations, the only 
versatile tool that can move across different scales and 
experimental embodiments. 

The generic theoretical system under study is depicted 
in Fig. 9: at normal incidence, a plane wave of frequency 
Lo excites a point dipole in the nanoparticle: 




p oc a„p(a;) 



(2) 



where the polarizability anp is described by the Mie 
theory [74, 75]. The poles of a„p define the optimal SERS 
frequencies. The re-radiated near-field from this dipole 
scales as r~'^, and is responsible for the enhanced absorp- 
tion. This will excite a Raman dipole; 



p' oc ani'^s-.'^y 



"-np 



{u) 



(3) 



where ajj(ws,w) is the Raman polarizability and uJs the 
Stokes-shifted emission frequency. This dipole near-field 
will in turn excite a secondary dipole in the nanoparticle: 



p" (X anp{ujs)r ^aR{uJs,uj)r ^a„p{uj) 



(4) 



now at the emission frequency lOs- Thus, the additional 
surface-enhanced Raman signal is: 



AIsERS'Xut / \p"\^dS 



(5) 



where the integration is over the SLG area. Assuming a 
square array of nanoparticles with spacing L and normal- 
izing to the corresponding Raman signal in the absence 
of the nanoparticles, we get the SERS enhancement: 



A/; 



SERS 



-aQicofQiu^sr (^) (6) 



/o 9 

where h is the separation between the nanoparticle cen- 
ter and the SLG plane, a = tto^/L'^ is the relative 
cross sectional area covered by the nanoparticles, and 
Q{uj) = |a„p(a;)| /47ra'^ is the the Mie enhancement. In 



FIG. 9: Scheme for SERS of a nanoparticle on graphene. 
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FIG. 10: FDTD simulations of suspended nanodisks (sym- 
bols) compared with with Eq. 6 (lines) for G peak enhance- 
ment. The nanoparticle radius a and elevation h in Eq. 6 
are adjusted to match the simulation peak enhancement. 
These are 2a=118nm, /i=32.2nm for the 140nm disks and 
2a=130nm, /i=38nm for the 210nm disks. Note that the to- 
tal enhancement factor F' is just F' = AIsers/Io + 1- 



general, Q depends on particle size and shape[75] (for 
spherical particles a„p(a;) is given by Eq. 8 in Meth- 
ods). In the limiting case of spherical particles with 
a <S^ X the size dependence in the Mie enhancement is 
removed and Q{io) = |[e(w) — l]/[e(w) + 2]|, where e{uj) 
is the nanoparticle's dielectric function. 

This is our main theoretical result: the Raman en- 
hancement scales with the metallic nanoparticle cross 
section, with the fourth power of the Mie enhancement, 
and inversely with the tenth power of the separation be- 
tween the nanoparticle centre and graphene. 

Since Eq. 6 does not take into account multiple reflec- 
tions in the substrate, we compare it with FDTD simula- 
tions of Au nanodisks suspended in air, shown in Fig. 10 
for the G peak (where there is a small blueshift due to 



the smaller effective refractive index below the nanoparti- 
cles without a substrate). In the absence of an analytical 
expression for the polarizability of a nanodisk, we use 
the polarizability of a sphere, shown in Eq. 8 in Meth- 
ods. To fit the simulation (symbols in Fig. 10) we adjust 
the radius a (which determines the peak wavelength) and 
separation h (which determines the peak value) in Eq. 6. 
As a result, the 140nm nanodisk is best represented by 
a sphere of diameter 2a=118nm elevated at ft,=32.2nm 
above the SLG plane, while for the 210nm nanodisk a 
sphere of diameter 2a=130nm at ft,=38nm is needed (lines 
in Fig. 10). These numbers are close to the actual values 
for the nanodisks. The agreement is good, considering 
the difference between disks and spheres, confirming that 
our theory captures the essential physics of SERS in 2d. 

Equation 6 is a valuable optimization/design tool. In 
particular, it identifies 3 steps to further improve SERS: 
1) larger nanoparticle coverage a, 2) larger Mie enhance- 
ment Q, 3) smaller nanoparticle-graphene separation h. 
The first is straightforward. The second is shape re- 
lated, e.g., ellipsoids have a different Q, which, for certain 
orientations, is stronger and red-shifted compared to a 
sphere[75]. This also influences the third step since, e.g., 
a flat oblate spheroid has a smaller distance h between 
its center and graphene. Thus, SERS enhancements can 
be much larger for ellipsoids and disks. This can already 
be seen in the comparisons made in Fig. 10. There, in 
order to match the disk simulation results, the represen- 
tative spheres had to be placed at an elevation h < a (i.e. 
so that they "cut" into the SGL plane). If we, however, 
re-evaluate Eq. 6 for the same representative spheres at 
h = a (i.e. so that they "rest" on the SLG plane), we lose 
more than two orders of magnitude enhancement. This 
is confirmed by detailed FDTD calculations. 

Finally, we note that the analytical expressions derived 
in Methods can be extended to flakes of increasing num- 
ber of layers, by vertical integration. Note, however, that 
if the number of layers is large reflection becomes impor- 
tant (it is 1% for 7 layers but exceeds 5% above 18 layers) 
and has to be taken into account. 



CONCLUSIONS 

We studied SERS in graphene patterned with a square 
array of Au nanodisks on SiO2(300nm)/Si. Significant 
enhancements were measured for both G and 2D bands 
at 633nm. Similar results were obtained for both disk 
sizes. Large-scale FDTD simulations reproduce well the 
experiments. To elucidate the physics of SERS in 2d, we 
derived analytic expressions, and showed that taking into 
account the SPR near-fields only, a simple closed-form ex- 
pression is found, where the Raman enhancement scales 
with the nanoparticle cross section, the fourth power of 
the Mie enhancement, and inversely with the tenth power 
of the separation between the nanoparticle centre and 



graphene. This points to thin nanodisks to achieve the 
highest SERS for 2d systems like graphene. 
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METHODS 

Let us consider a generic configuration as in Fig. 9, 
comprising a Au spherical nanoparticle of radius a at a 
distance h from SLG, normal plane wave incidence, with 
field amplitude Eq, frequency uj and polarization along 
X. The re-radiated fields from the nanoparticle are due 
to an induced electric dipole at its center: 



: eoa„p(w)E « eoanpioj)Ei 



'ox 



(7) 



where E is the local field at the nanoparticle, modified 
from the incident one due to the presence of the graphene 
layer. The SLG normal incidence reflectance is almost 
zero[71] and, due to continuity, the local fleld can be 
taken approximately the same as the incident one. This 
is further corroborated from calculations of the depolar- 
ization matrix of single wall nanotubes, where for po- 
larization along the nanotube axis, no depolarization is 
found[52, 76]. The nanoparticle polarizability a„p is de- 
scribed for a sphere in the Mie theory[74, 75]: 



,{uj) 



67ri h')pi{hka)')p'i{ka) ~ 'ipi(ka)')p'i{nka) 
k^ nipi(hka)£^[{ka) — ^i{ka)ip'iinka) 



(8) 



where ipi and ^i are the Riccati-Bessel functions, k — uj j c 
and h = h(iS) is the nanoparticle complex index of re- 
fraction. In general, the refractive index of metallic 
nanoparticles differs from the bulk value due the re- 
duced free electron relaxation caused by electron sur- 
face scattering[63, 77]. We have relatively large particles, 
where this correction is negligible. Thus we use the bulk 
Au refractive index. If the nanoparticle is much smaller 
than the wavelength, Eq. 8 simplifies [75]: 



a„p(w) w 47ra' 



;e(w) 



1 



e(cj) + 2 



(9) 



with e(w) = r^iiS) the nanoparticle dielectric function. 
The total field at at position r on SLG is [78]: 



Et(r,cj) «£;o(r,w)x- 



Akv 



47ren 



— f X (p X f ) 
r 



3- " l2 [3i"(r • P) - P] 



(10) 



The first term in the square bracket is the radiation- 
field, while the second is the near-field. They scale with 
distance and wavelength as (A^r)^^ and r^-^, (Ar^)^^ 
respectively. Assuming the nanoparticle distance from 
SLG to be much smaller than the wavelength, r <C A, 
and the resonant near-field much larger than the inci- 
dent field Eq, then the near- field r~^ has the dominant 
contribution. Its transverse component drives the SLG 
enhanced absorption. Simultaneously, the Raman emit- 
ted field also interacts with the Au particle SPR, further 
enhancing the total Raman emission. To fully explore 
both processes, we evaluate two quantities: 1) total ab- 
sorption enhancement; 2) total Raman enhancement. 



Absorption enhancement 

The additional absorption in graphene is due to 
the enhanced near field. Absorption is defined as 
current X field [78], thus approximated as: 



AA = Go \K \'dS 



P."/|2, 



(11) 



where Go — e^/AH is graphene's dynamical (optical) 
sheet conductance [71, 79-82]. Combining Eqs. 7,10, the 
near-field transverse component of the driving field is: 

|Ej;,^'H| = ii;ogMaV-3|[3f(f .i)-i]||| (12) 

with Q{uj) = \anp{Lu)\ /Aira^. In the above we ignore 
cross terms between incident and re-radiated fields. This 
is justified when strong near-field enhancements are ex- 
pected (as it should be for SERS), but less so when the 
near fields are of the same order as the incident field. 

From Fig. 9, a; = rsinO coscj), y — rsm9sm(j) and 
z — r cos 9, while r ~ h/ cos 6. The integration surface el- 
ement is pdpd(f) ~ h^ sinScos^'^ 6d(j)d6, where p — htanO 
is r's projection on the plane. Eq. 11 then becomes: 



2„6 



AA 



GqEqQ a 



2ir 



cos^ 9 sin ef{e,(t))d(t)de (13) 



with f{9,(j)) ==9 sin 



4^cos2 



AA 



— 6 sin^ 6 cos^ ( 
8h^ 



1. Then: 



(14) 



For nanoparticles arranged on a square lattice with spac- 
ing L, the absorption enhancement becomes: 
AA _ 37rQ2(^)a6 

where a = Tia^ / L"^ is the nanoparticle relative cross sec- 
tion. For spheres directly placed on SLG, i.e. h = a, 
Eq. 15 simplifies to: 



r^'^-Kl 



(15) 






(16) 



We remind here that in the a <C A limit, the Mie en- 
hancement is Q(w) = \[e{uj) — l]/[e(cj) + 2]|. 



Raman enhancement 

Going back to Eq. 12, E?! (w) will excite a dipole field 
on SLG at the Raman frequency ujs'. 



p' ^aii{u}s,u})\'Ell\p 



(17) 



The polarization of the Raman dipole is not necessarily 
the same as that of the driving field. Thus, for general- 
ity we assume this to be randomly polarized on the SLG 
plane. In this case it suffices to take the average of two 
dipoles, one polarized along x, and another along y. The 
one along x will emit as the dipole term of Eq. 10. We 
are again interested in the dominant near-field term that 
decays as r~^. This will get coupled to the nanopar- 
ticle and thus SPR enhanced. That is, it will excite a 
secondary dipole at the nanoparticle: 

p" = g(w,)aH(^.,w)Q(^)Soa''r"®|[3f(i-x)-x]|||2 (18) 

where we again consider the projection of the dipole given 
the backscattering geometry considered here. The radi- 
ated fiux can be taken as the additional surface-enhanced 
Raman signal: 



A/; 



SERS 



247reo 



p'"dS 



(19) 



where kg = ujg/c and we multiply by a factor 1/2 since 
we only consider the upper half fiux. Using the angular 
relationships for x, y, z and dS we get: 



AIsERS = 
r.27r 



cktEiQ^Lj)Q^{LJs)a^^ 
247reo/ii° 



aR{ujs,uj)\'^ X 



where fa;{0, 

4 



cos^e sine f^{e,(f))d(f)de (20) 

"'0 

<j)) = 81sin®0cos4 - 108 sin*^ 6* cos" + 
0(1 + 2cos2 0) - 12sin2 6'cos2 + 1. The 
angular integration yields 337r/280. For the calculation 
with the Raman dipole along y, the angular part of 
Eq. 18 becomes |[3f(f • x) — x]||||[3f(f • y) — y]||| with a 
slightly different /y(^, ci), but a similar angular integra- 
tion value of 277r/28G. The average angular contribution 
is 37r/28 « tt/9. Thus the additional Raman signal can 
thus be written as: 

^^^^«^"9 24W^I5 l"«(-.,-)l (21) 

To evaluate the enhancement factor, we normalize to 
the expected signal /o in the absence of the nanoparti- 
cles. Considering a square unit cell of side equal to the 
nanoparticle spacing L, we get: 



Jo«L^ 



, c/c! R 



-\aR{uJs,uj)\ 



247reo 
The Raman enhancement factor is 

AIsERS ^ 

h ^ 9 



(22) 



-<yQ\u)Q\us) (-) (23) 



The Finite-Difference Time-Domain Method 

In the FDTD method, MaxweU's equations are time- 
integrated on a computational grid: 



V X E = -^idtH 



(24) 



N 



VxH = eoe^dtE + dtPo + J2dtPj (25) 

where material polarization is taken into account through 
the polarizabilities P: 

d^Po+-/dtPo = WpfoE (26) 

dfPj+TjdtPj+npj = Ae^fl^eoE (27) 

This gives a Drude-Lorentz model for the dielectric 
function [83]: 



eic) 



w; 



w^ + iuj"f 



N 

3 = 1 "j 



Ae,^,^ 



j'^ — iujTj 



(28) 



where the first term is the Drude free-electron con- 
tribution and the second contains Lorentz oscillators 



corresponding to interband transitions. 



and 1/7 



are the free electron plasma frequency and relaxation 
time, flj, Aej, r,„ are transition frequency, oscillator 
strength and decay rate for the Lorentz terms. To ac- 
curately reproduce the experimental dielectric functions 
(Au and Cr from Ref.[68], SLG from Rcf.[69]) we treat 
these as fit parameters. For Au we use N=A, and 
£^=3.454, Aej=(0.376, 0.63, 1.208, 1.124), fia;p=8.73eV, 
n7=0.046eV, n%=(2.72, 3.13, 3.88, 4.95)eV, and 
hTj={0.39, 0.655, 1.16, 1.67)eV. For Cr we use A^=3, and 
£00=1, Aej=(9.54, 15.5, 1.1), fia;p=5.51eV, n7=0.731eV, 
hnj={lA3, 2.36, 3.64)eV, and hTj={l.l9, 1.94, 1.41)eV. 



Finally, for SLG we use N=3, and eoo=l-964, Ae, 



.99, 



1.69, 1.53), ;za;p=6.02eV, n7=4.52eV, nrjj=(3.14, 4.03, 
4.59)eV, and hrj={7.99, 2.01, 0.88)eV. Fig. 11 plots our 
model dielectric functions along with the experimental 
ones, showing an excellent agreement. 
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